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On the Relation between Period and Density of Algol-Variables. 
By the Rev. J. Stein, S.J., Sc.D. 

(1 Communicated by Prof. H. H. Turner t D. Sc., F.R.S.) 


1. As is well known, a maximum value of the mean density 
of an Algol-system can b$ derived from the period (P) and the 
total duration of eclipse (2t 0 ).* If the orbit is supposed to be 
circular, this maximum-value I) is given by. 


D 


K ■ n = — 
P 2 sin s nt 0 ’ F ’ 



where P and t 0 may be expressed in hours, D in the mean density 
of the Sun as unity. 

This value is identical with the real mean density (S) if the 
two stars are of the same size, and if the inclination of the line of 
sight to the orbit is zero. 

In order to determine the constant K we put— 


P = one year = 365*25 x 24 11 ; 7^ = 32' 3"*64 = mean apparent 
diameter of the Sun; consequently D = £ and K = 31*17. 

Thus 


D = 3 11 ? 

P 2 sin 3 nt 0 


(*> 


and 31*17 being =1*005 7r 3 , we can bring this into the simple 
form— 


D 


P / nt 0 V 

(2 1 0 ) 3 \sin ntj 


( 3 ) 


M6riau has shown that D is not very different from S, if one 
star is not considerably larger in size than the other. 

2. In the Mitteilungen der Hamburger Stemwarte , No. 11, 
Dr. Graff has deduced from his own observations the elements of 
the orbits of 10 Algol-variables. A slight extension was given to 
these by Professor Ristenpart ( Ast . Nach ., No. 4250), who derived 
from the elements the mean density of the systems by the formula 


s = a3 ° 

P 2 (i+k 3 ) 

where a is the radius of the relative orbit, k the radius of the dark 
satellite, both expressed in the radius of the bright star as unity. 
If P is given in hours, then C = JK. 

Arranging the stars according to decreasing periods, Ristenpart 
finds a nearly progressive increase of density; and he adds that 
this might be expected with regard to formula (4): “ NatiLrlich 

* M. M6riau, “ Density des etoiles variables du type d*Algol,” Comptes 
Rendus de VAcad. d. Sciences, vol. 122 (1896), p. 1254. 
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laufen die Korper um so rascher um einander um, je dichter sie 
sind.” 

This argument is, of course, understood, “cseteris paribus”; 
strictly we can only infer from the results obtained that, with 

eft 

increasing P, the factor -- in general increases less rapidly 

i + K 3 


than P 2 . 

3. In order to examine more closely the relation between density 
and period, I have put together the following list of 38 Algol-stars, 
for which I could find the duration 2 1 0 of the eclipse, determined 
exactly enough for the purpose. The fourth column contains the 
maximum-density D; the next one gives the mean densities 8, 
derived by Ristenpart from Graff’s elements. The slight differences 
D - 8 are an experimental proof that D is a fairly approximate 


value of 8. In the sixth column are 
( = “ rarity ”). 


given the values of = R 


Star. 

P. 

2 t 0 . 

D. 

8. 

K,. 


h 

h 




RZ Draconis (Z 2 ) 

13-2 

27 

0*848 


1*2 

ZZ Cygni (Z 6 ) 

15*1 

3*5 

•358 


2*8 

U Ophiuchi 

20*1 

5 *i 

*211 


4*8 

RX Herculis (X 2 ) 

21 *4 

5 

•226 


4*4 

RU Monocerotis (U 2 ) 

21 *6 

5 

•228 


4*4 

R Canis Mai. 

27-3 

5 

•259 


3*9 

RZ Cassiopeiae (Z 2 ) 

287 

5*5 

•203 


4*9 

Z Draconis 

32 *6 

47 

*347 

*345 

2*9 

RX Draconis (X 2 ) 

45*5 

5 

•389 


2*6 

RW Monocerotis (W 2 ) 

45*8 

6*5 

•190 


5*3 

RV Persei (Y 2 ) 

47*4 

6*5 

•194 


5*2 

5 Librae 

op 

vn 

ui 

12 

O 4 I 


24*4 

RS Sagittarii (S 2 ) 

58*0 

10*4 

*061 


16*4 

U Cephei 

59*8 

10 

•069 


, 14-5 

RY Aurigae (Y a ) 

65*4 

8 

•138 


7*3 

143*1907 Andromedae 

66*4 

8 

•140 


7 * 1 

RW Tauri (W 3 ) 

. 66*4 

7*9 

•148 

•136 

6*8 

RW Geminorum (W 2 ) 

68*8 

12 

*046 


23*9 

£ Persei 

68*8 

10 

*0 77 


13*0 

Y Cygni 

71-9 

8 

•150 


67 

Z Persei 

73*3 

li-i 

•060 

•059 

167 

Y Camelopardalis, 

79*3 

12 

•051 


197 

' WW Cygni (Z 5 ) 

79*6 

n*8 

•053 

•052 

18 *9 

U Sagittae 

81 r i 

13-1 

•041 

*041 

24*4 

UW Cygni (Z 4 ) - ’ 

82*8. 

io*5 

•079 

*075 ' 

t 127 
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Star. 



D. 

8 . 

It. 


h 

h 




U Coronse 

82*8 

97 

•097 


10*4 

A Tauri 

94*9 

10 

*IOO 


10*0 

Z Herculis 

95*8 

9*2 

*129 


7*8 

R Arte 

106*2 

10*3 

•103 


9*8 

SW Cygni (V 8 ) 

109-5 

ii*8 

•070 

*064 

14*3 

RR Delphini (R 2 ) 

1104 

14 

•043 


23*3 

W Delphini 

ii 5*3 

17*2 

•025 

•025 

40*0 

S Velorum 

142*4 

15*2 

*043 


23*3 

SY Cygni (X 3 ) 

, i 44 *i 

19*0 

•023 

•023 

43*5 

RY Persei (Y 2 ) 

164*7 

26 

•010 


100*0 

YW Cygni (Y 5 ) 

2 02 3 

20 

■027 

*020 

37 -i 

S Cancri 

227*6 

’21*5 

•025 • 


40*0 

RR Puppis (R 2 )' 

250*3 

1 7 

*033 


3°‘4 


4. A rough preliminary calculation showed that there can be 
no question of proportionality between P 2 and R. On the other 
hand, the product P D oscillated between relatively narrow limits. 
This seemed to make an investigation of the correlation between 
P and R more promising. 

By taking direct square and product deviations from mean 
values we find :—* 

Including all the 38 stars— 


Mean period P =82*66 hours 

Variability’in period <r P = 55*60 hours 
Coefficient of variability V P = *672 


Mean *‘rarity ” B =1695 
Variability in rarity o* R = 16*28 
Coefficient of variability V R = *960 


Correlation : 


nt,p='733±'°5 I ( P - e •)• 


Probable value of rarity : 


R= -215 P - -83 ± 7-46 (5) 


Excluding 4 stars with period greater than 6 days, we have for 
34 stars— 


Mean period P =61*65 hours 

Variability in period <r P = 34*11 hours 
Coefficient of variability V P = *554 


Mean “ rarity ” R =12*86 
Variability in rarity or R = 1 o * 13 
Coefficient of variability V R = *788 


Correlation: 

?r,p = *73i ±*°54 


Probable value of rarity : 

R= *217 P - *51 ± 4.66. 

* For the notations used, see M.N ., Dec. 1908, vol. lxix. pp. 128 ff.. 
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5. Thus the magnitude of correlation seems to indicate a real 
relationship between period and density, or, more exactly, between 
period and duration of eclipse (see Note at end). For a longer 
period (say > 3 days) we can then put with sufficient accuracy: 


or, according to (2): 


PD = = 4-63 

210 


P sin 3 nt Q — 673. 


When P is very large, nt 0 must be small, hence nt 0 = sin nt 0 , and 


(2jy 3 =-2i6 P 2 . . . . . (6) 

Applying this, for instance, to the Algol-star 142, 1907 Cassiopeise, 
with a period of 36^56 = 877*4 hours, we find : 


2 ^o“55 h° urs i D = 0*005. 

And, in fact, according to Euebo ( Ast. Nach., No. 4241), the dura¬ 
tion of eclipse lies between 49*7 and 64*8 hours (0*003 < <0*007). 

It is hoped that this result, in connection with the relations 
found by Professor Karl Pearson and Miss Julia Bell, may throw 
some new light on the nature and evolution of this most interesting 
class of variable stars. 


Specola Vaticana: 
1909 January 28. 


Note ,—As was remarked in No. 5, the relation found between 
P and D is, properly speaking, a relation between the period and 
the observed duration of eclipse. This point was especially 
emphasized by Professor Karl Pearson, to whom the MS. was 
submitted by Professor Turner. He makes the following interest¬ 
ing remark:— 

“With reference to the formula 

D= 3 i ‘*7 . „ = 

P 2 sin 3 (nt 0 ) ’ P’ 

if t 0 be a given function of P, it follows that D would be absolutely 
determined by P. The fundamental question, therefore, is what is 
the relation between period and eclipse time ? I think these must 
be very closely related. I have accordingly worked them out, and 
find 

rp,t n = +‘877. 

“It will thus, I think, be seen that D being a mathematical 
function of P and t 0 , not the true density, the correlation of D 
and P is a secondary result of the fundamental relation of the 
observed quantities P and £ 0 . In fact, the value of r EP can be 
deduced at once from the values of P, i 0 , crp, <rt 0 and 
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According to Professor Pearson’s remark, putting together the 
elements of correlation between P and 1 0 , we get— 


Mean period P =82*66 

Variability in period cr P = 55*60 
Coefficient of variab. V P = *672 

Correlation : 


Mean Eclipse time 2t 0 — 10*532 
Variab. in m.E.t. o- 2 * 0 = 5*236 
Coeff. of variab. V 2 f 0 = *497 


rp,2«o = "877 ± '025 (p.e.). 
Probable value of eclipse time, 

2t 0 = *0826 P + 3*70± 1*70 (p.e.) 


(7) 


This value of 2* 0 , between the limits considered, will, of course, 
not differ substantially from that given by the formula (5): 


-o P 2 sin 3 nt n ^ 

R-=0*21 5P - 0*83 . 

31*17 D . ** 


• • (5) 

This is easily shown by means of the approximate formula (6'): 

2 ^o = °’594 P 2/3 .(6') 

In fact, if we wish to replace the curve 
‘ y = Ci» 2 /s; a < x < p 

approximately by a straight line 

y = a + bx; a <x < /3, 

the best fitting values of a and b , according to the method of least 
squares, are given by the equations 

a^dx 4 * b^xdx = c^^dx 

ajxdx + x 2 dx ~ cjxVsdx 

where the integrals are to be taken between the limits a and /?. 
Putting y = 2 1 0 , x = P, 0=0*594, a =25. /3= 250, we find 

2^ = °*°797 P + 4 * 5 ° 
in close agreement with (7). 

The following table gives in the 2nd. and 3rd. columns the values 
of 2derived from formulse (5) and (7): 


p 

(2*0)5. 

(2*0)7- 

Diff. 

h 

h 

h 

h 

2 5 

5 ‘ 2 

5*8 

- *6 

5 ° 

8*5 

7*8 

+ 7 

75 

11*0 

9*9 

+1 * I 

IOO 

13*4 

12*0 

+ i '4 

* 5 ° 

17*4 

16*1 

+ i *3 

200 

20*8 

20*2 

+ *6 

250 

24*4 

24*3 

+ *i 
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Radial Movement in Sun-spots. By J. Evershed. 

Displacements of the absorption-lines of hydrogen and calcium 
are frequently observed in all spot disturbances where there is 
active change in progress, and during the genesis of a spot. 
These line-shifts may amount to several Angstrom units in either 
direction in the spectrum, indicating velocities towards or away 
from the observer of one or more hundred kilometres per second. 
Only the elements of the higher chromosphere seem to be affected 
as a rule, and the movements are seldom maintained for more than 
a few minutes at a time; also they are rarely found within the 
penumbral area of a spot. 

Recently line-shifts of quite another character have been 
photographed here. These are permanent, affect most, if not all, 
of the Fraunhofer lines, and are greatest in the penumbrse. 

In some spectra obtained by me this year, in the fourth and 
fifth orders of a Rowland plane grating of 3*2 inches ruled surface 
and 14,428 lines to the inch, the lines crossing the spot-band 
were found to be slightly bent, or inclined one or two degrees as 
compared with the lines in the neighbouring photosphere. The 
first plate showing this feature clearly was exposed on January 
7, 1909, in the fourth order of the grating. It covers the region 
A4650 to A4790, and the scale is iA=i*o 8 mm. The total dis¬ 
placement of the lines measured on each side of the penumbra 
amounts to 0*02 7A, indicating a receding velocity on the north¬ 
west side of the spot of 0*85 kra. per second, and an approaching 
velocity of the same amount on the south-east side. The spot was 
in latitude 9 0 north and longitude 31 0 west of the central meridian, 
and the position-angle of the north end of the spectrograph slit 
was 314°. 

The appearance of the lines on this plate at once suggested a 
rotation of the absorbing gases in the penumbra—not a vortex 
motion, but a rotation of the spot as a whole about a point at 
its centre. The lines are quite straight over the spot, although 
inclined to the undisturbed lines; and in the centre of the umbra 
there is no displacement. The hypothesis of circular motion of 
any kind has proved, however, to be certainly untenable. From 
an examination of about 150 spectra, obtained since January 7, 
and representing seven different spots in the northern hemisphere 
and four in the southern, the following definite results have been 
obtained:— 

1. All the spots examined show line-shifts of about the same 
order of magnitude when at the same distance from the centre of 
the Sun’s disc. 
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